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ABSTRACT 

Based on analysis of photometric observations of nearby M type stars obtained with ASAS, 31 periodic 
variables were detected. The determined periods are assumed to be related to rotation periods of the inves- 
' I ' tigated stars. Among them 10 new variables with periods longer than 10 days were found, which brings the 

— total number of slowly rotating M stars with known rotation periods to 12 objects. 

X-ray activity and rotation evolution of M stars follows the trends observed in G-K type stars. Rapidly 
Q 1 rotating stars are very active and activity decreases with increasing rotation period but the period-activity 

' relation is mass-dependent which suggests that the rotation period alone is not a proper measure of activity. 

^ ' The investigated stars were grouped according to their mass and the empirical turnover time was determined 

' for each group. It increases with decreasing mass more steeply than for K type stars for which a fiat de- 

pendence had been found. The resulting Rossby number-activity relation shows an exponential decrease of 
activity with increasing Rossby number 

The analysis of space motions of 27 single stars showed that all rapidly rotating and a few slowly rotating 
^ ' stars belong to young disk (YD) whereas all old disk (OD) stars are slowly rotating. The median rotation 

' period of YD stars is about 2 days and that of OD stars is equal to 47 days, i.e., nearly 25 times longer. The 

^ average X-ray flux of OD stars is about 1.7 dex lower than YD stars in a good agreement with the derived 

ly-^ , Rossby number-activity formula supplemented with rotation-age relation and in a fair agreement with recent 

■ observations but in a disagreement with the Skumanich formula supplemented with the activity-rotation 



relation. 
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'. 1 Introduction 

k>( \ Chromospheric-coronal activity is inherent to lower main sequence (MS) stars. In ro- 

}_j ' tating stars with subphotospheric convection zones magnetohydrodynamic dynamo op- 

I erates generating magnetic fields inside the zone and/or at the interface between the 

convection layer and the radiative core. After emerging above the photosphere, the 
magnetic field drives the heating of a stellar atmosphere and produces various activity 
phenomena like spots, plages, flares and magnetized winds. Empirical data indicate 
that the activity level increases with increasing rotation rate (Hartmann and Noyes 
1987, Maggio et al. 1987, Stgpien 1989, 1994, Hempelmann et al. 1995). The ob- 
served period-activity relation for MS stars from a narrow spectral type shows a small 
scatter but when stars of all spectral types are plotted together, the scatter increases sub- 
stantially (Noyes et al. 1984, Stgpieri 1994, Pizzolato et al. 2003). This clearly shows 
that period-activity relation varies with stellar mass. To decrease the scatter, rotation 
periods of stars of a given spectral type must be scaled down by a mass dependent 
quantity called the convective turnover time Xc- Such a scaling is suggested by a simple 
parametric dynamo theory (Durney and Latour 1978, Gray 1982) in which the Rossby 
number Ro = Pmt/ rather than rotation period is a primary parameter controlling the 
efficiency of field generation. 

Values of turnover times for stars of various masses can be determined from the- 
oretical models (e.g.. Oilman 1980, Oilliland 1985, Rucinski and VandenBerg 1986, 
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Kim and Demarque 1996) or empirically (Noyes et al. 1984, St§pieri 1994, 2003, Piz- 
zolato et al. 2003). The values of Xc for G-K type stars, given by various authors, 
do not differ much but for M type stars they diverge badly. Theoretical models indi- 
cate a steep increase of with decreasing mass. For example, the theoretical value of 
turnover time for a 0.5 M© star is 3-3.5 times longer than for a 1 M© star (Kim and De- 
marque 1996, Rucinski and VandenBerg 1986 - in the latter case a sUght extrapolation 
is needed because their models reach only 0.7 Ma,). Purely empirical determinations 
indicate that Xc increases with decreasing mass down to about 0.8 Mq, but then it levels 
off till about 0.6 M© (St§pieii 1994, 2003, Pizzolato et al. 2003). Due to a shortage 
of known periods for less massive, slowly rotating stars, the behavior of the empirical 
turnover time with decreasing mass beyond 0.6 M© is very uncertain. There are indi- 
cations that it may start increasing again but the situation has been obscure. The results 
of a recent analysis of angular momentum evolution of very low mass stars (VLMS), 
based on observations of young clusters, suggest that very long turnover times, pre- 
dicted by theoretical models, do not agree with the observed rotation velocities of such 
stars. Sills, Pinsonneault and Terndrup (2000) noted that the relation (£)crit{Xc), where 
cOcrit is a critical angular velocity for the occurrence of saturation, breaks down if val- 
ues of Xc for low mass stars are extrapolated from the higher mass models of Kim and 
Demarque (1996). Lower values are in a better agreement with observations. 

M type stars seem to share all characteristic properties of G-K type active stars: the 
most active stars have Lx/Lboi at the level of 10~^ i.e., close to the upper limit for all 
active stars whereas the lowest measured X-ray fluxes are some 3 orders of magnitude 
lower as can be seen in the NEXXUS database (Schmitt and Liefke 2004). High activity 
level is related to stellar youth and it decreases with age so that old stars show notably 
decreased level of activity, both chromospheric (Delfosse et al. 1998, Silvestri, Hawley 
and Oswalt 2005) and coronal (Fleming, Schmitt and Giampapa 1995, Feigelson et al. 
2004). The most active M stars rotate rapidly, with periods from a fraction of a day up 
to several days. Quite a number of rotation periods of such stars is known (Pizzolato et 
al. 2003 and references therein). Low activity stars seem to rotate slowly - they show 
values of vsini below the present resolution threshold of about 2-3 km/s (Delfosse et 
al. 1998), which corresponds to a lower limit of several days for the rotation period. 
This is not a very restrictive condition because the lowest mass stars with such rotation 
periods are still in a saturation regime (Pizzolato et al. 2003). To extend the period- 
activity relation to the least active M type stars we need direct measurements of rotation 
periods of slowly rotating stars with periods several times longer than the saturation 
Umit. 

A standard method of measuring rotational periods is based on observations of 
stellar surface inhomogeneities which produce signal modulations over the rotational 
period. Several periods of chromosphericaUy active stars have been determined in this 
way from rotational variations of Call emission, carried out within the Mount Wilson 
program (Noyes et al. 1984, Donahue, Saar and Baliunas 1996 and references therein). 
The technique works well for F-K type stars. In M type stars the emission in Call lines 
becomes a poor measure of the chromospheric energy losses and is difficult to measure. 
So far, the only M type star for which a rotation period has been determined from 
calcium emission modulation, is GJ 41 1 (=HD 95735) with spectral type M2V (Noyes 
et al. 1984). Star-spots produce another type of surface inhomogeneities, observable 
with a broadband photometry. A rotational modulation of stellar brightness has been 
observed in many active stars, including those of M type (see Messina et al. 2003 for 
a recent review). High activity stars have typical amplitudes of several percent which 
is easy to detect, particularly for short rotation periods. Low activity stars are expected 
to have much less spotedness producing correspondingly weaker light variations and 
on a much longer time scale (Messina et al. 2003). In addition, possible variations of 
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a spot pattern on a time scale comparable to the rotation period can mask rotational 
modulation. Long time series of observations are needed to filter out the correct value 
of a period. Because of all these problems, only very few detections of rotation periods 
of slowly rotating stars have been reported from the observations of star-spots. 

In the last decade many very low mass stars (VLMS) and brown dwarfs have been 
detected using high sensitivity observational techniques in red and infrared. Several 
studies followed, including searches for rotation periods and measurements of vsini. 
Most of the studies are concentrated on young clusters (Scholz and Eisloffel 2005 and 
references therein) but several field objects have also been observed. The results show 
that many observed VLMS are active and fast rotating. The transition from partly 
to fully convective stars has no apparent influence on their activity and/or rotation 
(Delfosse etal. 1998, Mokler and Stelzer 2002). The average value of vsini apparently 
increases with the advancing spectral type from early M to late M and L (Delfosse et al. 
1998, Mohanty and Basri 2003, Scholz and Eisloffel 2005). This is usually interpreted 
in terms of increasing time scale for spin down with decreasing mass but a detailed 
age-activity-rotation relation for M type stars of various masses is unknown. One 
should not forget, however, when interpreting observations of stars with ages of about 

1 Gyr or less, that evolutionary effects should also be taken into account due to a rapid 
increase of the approach time to ZAMS of VLMS, up to 3 Gyr for least massive stars 
(Baraffe et al. 1998). More stars with known age, activity level and rotation period 
(particularly those inactive and/or slowly rotating) will help in solving the problem of 
spin down and activity decrease of low mass stars. 

The present investigation is aimed at enlargement of the number of M type stars 
with known rotation periods. A particular care is taken to detect slowly rotating stars. 
To this purpose an extensive data set of photometric observations obtained during 
the automated sky survey ASAS (Pojmanski 1997, 2004) is used. To assure a min- 
imum reasonable accuracy, we analyzed M stars brighter than a visual magnitude 
y = 12.5 mag which restricted our sample to M0-M3.5 type stars only, with an ex- 
ception of Proxima (M5.5). 

Section 2 describes the star selection criteria and method of analysis. In Section 3 
the results of the period search are presented and interpreted in terms of empirical 
turnover time and age-activity-rotation relations. Section 4 contains the summary of 
the main conclusions. 

2 Star Selection and Analysis 

ASAS cameras are located at Las Campanas Observatory, Chile, and the observations 
cover the sky south of declination -h28°. At the angular resolution of 14"/pixel several 
visual binaries are not separated. 

The most accurate ASAS photometric measurements are obtained for stars with V 
magnitudes between 8.0 and 12.5. Brighter stars are usually saturated, whereas the 
observational errors increase rapidly for fainter stars. For the period search we selected 
stars with the following properties: 

(i) X-ray flux is known from Hiinsch et al. (1999) or from NEXXUS database 

(ii) spectral type is M, 

(iii) 8 < y < 12.5, where y is a mean V magnitude. 

180 stars from the ASAS database fulfilled these conditions. Only good quality 
photometric data (with quality grade "A" and "B" - Pojmanski 2006, private commu- 
nication) were included into the period analysis. 

Objects monitored by ASAS are typically observed once per night or, sometimes, 
even more sparsely, depending on weather. The Nyiquist frequency corresponds in this 
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case to the period of 2 days. We looked for periods from 0.2 d up to the length of an 
observational run in a single season (w 100 days). We assumed that season to season 
variations are not connected with rotation and we eliminated them before performing 
period search. Because the shortest periods, found in the search, are shorter than the 
Nyiquist limit, it is probable that some of them may be aliases of the true periods. 

Many late type stars flare from time to time with typical duration of a flare of less 
than one day. Such events should be eliminated before period search, because points 
with large deviations from the mean value have large weights in period searching algo- 
rithms. Unfortunately, the ASAS observations are not well suited for identification of 
stellar flares because they usually manifest themselves as single data points indicating 
a substantial brightening. To get rid of them we rejected data points deviating from 
the mean seasonal magnitude by more than 3.5 standard deviations, or from the grand 
mean value by more than 0.35 mag. 

Apart from sudden brightenings many similar light drops occur within the data. It 
is not clear what causes them. The simplest explanation relates them to wrong mea- 
surements. Magnitude measurements fainter than the mean value are more frequent 
than those above the mean (Pojmanski 2006, private communication). We applied the 
same procedure to them as in case of brightenings. Altogether, we rejected 229 drops 
and 87 brightenings, which is about 0.6% of all data points. Only one star, GJ 551 (= 
Proxima) happened to brighten over two consecutive observations: first by 1.13 mag 
above the mean level and 1 .5 hour later it was still brighter by 0.35 mag. In other cases 
we do not have neighboring observations indicating significant brightness deviations 
from a normal level. There were on average 285 measurements per star left, extended 
over 6 seasons. 

Late type active stars often show season-to-season light variations with amplitudes 
of the order of 0.1 mag. A good example of such star is AB Dor which was observed 
nearly continuously over more than 20 seasons (Kiirster et al. 1997, Jarvinen et al. 
2005). The long term variations are usually interpreted as resulting from activity cy- 
cles. Several stars from our sample also show seasonal variations. Because of that 
we have done period search separately for each season unless the number of useful 
observations in the particular season was lower than 50 data points. Data from such 
seasons were merged with the ones from neighboring seasons and analyzed together. 
A similar period search was also performed on the whole data set after eliminating 
season to season variations. We used the AoV algorithm, together with the subrou- 
tine CAOV which can be downloaded from the Alex Schwarzenberg-Czerny homepage 
\http://w ww. camk. edu.pl/^alex^ A detailed description of the program can be found in 
Schwarzenberg-Czerny (1989). In short, we are folding and binning data points with a 
trial period, and for a given number of bins and data points we obtain AoV statistics 

®Aov^sysi (1) 

where 

{r-l)s\^Y^ni{Ti^x)^ (2) 



!=1 



and 



(«-'->2=EEN--^)' (3) 



with r being a number of bins, n,- - number of data points in the "/-th" bin, xj - mean 
value of data points in the "/-th" bin and x- mean value of all data points. We used six 
phase bins and assumed that the detected periodicity is significant when ©aoV (P) > 8.0, 
unless the indicated period was very close to 1 day due to a slow, systematic drift of the 
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stellar brightness over the whole season. Such periods were treated as spurious. The 
results of the period search are presented and discussed in Section 3. 

3 Results and Discussion 
3.1 Period Search 

Out of 180 analyzed stars only 31 showed statistically significant periodic variations. 
Table 1 lists basic data for them and the results of the period search. Spectral types, 
absolute magnitudes and X-ray data are taken from NEXXUS. For X-ray variable stars 
average values of logRx are taken. Bolometric corrections were calculated from color 
indices R-I, using a formula from Delfosse et al. (1998): BC = -0.083 - 0.501 (/? - 
7) -0.646(7? -7). 

Table 1 

Basic data for stars with significant periodic variability 



name 


ASAS 


d 


Mv 


Sp 


(logLi,) 


(logR.) 


P 


©AoV 


amp 


mass 




des. 


pc 


mag 








day.s 




mmag 


Mo 


GJ 84 


020504-1736.9 


9.4 


10.32 


M2.5 


32.08 


-4.4 


44.51 


24.4 


26 


0.45 


GJ 103 


023422-4347.8 


11.5 


8.55 


MO 


32.53 


-2.8 


1.563 


64.7 


63 


0.65 


GJ 1054A 


030755-2813.2 


18.0 


8.96 


MO 


32.36 


-3.0 


0.513 


15.4 


18 


0.63 


HIP 17695* 


034723-0158.3 


16.3 


10.53 


M3 


31.98 


-3.0 


3.880 


25.1 


39 


0.43 


Gl 176 


044256+1857.5 


9.4 


10.11 


M2 


32.10 


-4.7 


38.92 


19.5 


34 


0.48 


GJ 2036A 


045331-5551.6 


11.2 


10.89 


M2 


31.95 


-3.0 


0.849 


9.4 


11 


0.38 


GJ 182 


045934-^0147.0 


26.7 


7.96 


M0.5 


32.75 


-3.1 


4.410 


19.9 


27 


0.69 


GJ 3331A 


050650-2135.1 


11.8 


9.93 


M1.5 


32.10 


-3.0 


0.341 


8.4 


9 


0.50 


GJ 205 


053127-0340.6 


5.7 


9.18 


M1.5 


32.37 


-4.7 


33.61 


10.3 


8 


0.60 


GJ 3367 


054717-0000.8 


24.4 


9.05 


MO 


32.31 


-3.4 


12.05 


8.9 


14 


0.62 


GJ358 


093946-4104.1 


9.5 


10.86 


M2 


31.94 


-3.9 


25.26 


21.3 


14 


0.39 


GJ375 


095834-4625.5 


15.9 


10.26 


M3.5 


32.30 


-3.1 


1.877 


11.0 


22 


0.37 


GJ382 


101217-0344.7 


7.8 


9.81 


MI.5 


32.18 


-4.7 


21.56 


21.6 


11 


0.52 


GJ431 


113146-4102.8 


10.5 


11.42 


M3.5 


31.81 


-3.5 


14.31 


10.5 


31 


0.33 


TWA5A 


113156-3436.5 


50 


7.88 


M1.5 


33.04 


-3.1 


0.7767 


21.2 


36 


0.70 


GJ494 


130046-H222.5 


11.4 


9.46 


M0.5 


32.25 


-3.4 


2.889 


15.4 


16 


0.56 


GJ 551 


142942-6240.8 


1.3 


15.49 


M5.5 


30.76 


-3.8 


82.53 


24.0 


21 


0.11 


GJ 569A 


145429-H 606.1 


9.8 


10.24 


M3 


32.05 


-3.7 


13.68 


8.3 


10 


0.46 


GJ 9520 


152152+2058.7 


11.4 


9.83 


M1.5 


32.19 


-3.3 


0.369 


8.4 


13 


0.51 


GJ 618A 


162003-3731.7 


8.5 


10.96 


M3 


31.89 


-4.7 


56.52 


9.5 


11 


0.38 


GJ 2 123 A* 


165648-3905.7 


14.6 


10.36 


M3 


32.03 


-2.5 


0.320 


65.9 


46 


0.45 


GJ 669A 


171954+2630.1 


10.7 


11.28 


M3.5 


31.87 


-3.2 


0.950 


10.3 


30 


0.35 


GJ674 


172839-4653.7 


4.5 


11.10 


M3 


31.80 


-4.2 


33.29 


9.5 


8 


0.37 


GJ729 


184949-2350.2 


3.0 


13.09 


M3.5 


31.21 


-3.5 


2.869 


15.3 


10 


0.20 


GJ 799 A/B 


204151-3226.1 


10.2 


10.94 


M4.5 


32.17 


-3.0 


0.7813 


16.8 


13 


0.38 


GJ 803 


204509-3120.5 


9.9 


8.82 


Ml 


32.52 


-2.9 


4.848 


30.5 


46 


0.64 


GJ 1264A 


214905-7206.1 


16.1 


8.76 


M0.5 


32.53 


-3.3 


6.669 


68.5 


21 


0.64 


Gl 841A 


215741-5100.4 


16.2 


9.31 


M2.5 


32.50 


-3.2 


1.124 


12.4 


10 


0.58 


GJ 867A 


223845-2037.3 


8.6 


9.42 


Ml. 5 


32.36 


-3.2 


4.233 


11.1 


10 


0.57 


GJ890 


230819-1524.6 


21.8 


9.18 


MO 


32.24 


-3.0 


0.431 


11.0 


37 


0.60 


GJ 897A 


233247-1645.4 


15.4 


10.01 


M2 


32.22 


-3.1 


4.828 


12.0 


14 


0.49 


GJ411 


110321+3558.2 


2.54 


10.46 


M2 


31.89 


-4.9 


48.0 






0.44 


GJ699 


175749+0441.6 


1.8 


13.25 


M4 


31.17 


-5.4 


130. 






0.20 



Consecutive columns give names, ASAS designations (= coordinates given in the order: RA-Dec), distances, absolute 
visual magnitudes, spectral types, bolometric luminosities, logR^ = log(L,/Lb<)i), detected rotation periods, values of 
the parameter ®aov at the Usted period, amplitudes of light variations, and estimated stellar masses. 



For two stars with no 7? — 7 indices marked with asterisk the bolometric corrections 
were calculated from the absolute visual magnitudes, as suggested by Pettersen (1983): 
SC=-0.397Mv + 2.386. 

Stellar masses were calculated from the absolute magnitude My using the formula 
derived by Delfosse et al. (2000): log(M/MQ) = 10-3(0.3 + l.SlMy + 7.614M^ - 
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1.698M^ + 0.060958M4). 

In case of tight binaries, the magnitudes given in Table 1 were not corrected for 
fainter components, so their masses may be overestimated. The exception is GJ 375 
which consists of nearly identical components (Montes et al. 2006). The formula given 
by Delfosse et al. (2000) applies only to stars with My > 9 mag. For a few brighter 
stars the masses were found using the caUbration given by Andersen (1991). The last 
two entries in Table 1 are taken from literature. 

Left panels of Figs. 1-4 present all photometric measurements of periodic stars 
vs. hjd = HID — 2450000 d with observations showing periodic variations marked as 
filled symbols. Right panels show filled symbols folded with the detected periods after 
applying corrections for seasonal variations. In most cases periodic variability is seen 
over the whole observational interval without any significant phase shift from one sea- 
son to another, but a few stars show transient periodic variability. The number of filled 
symbols and the total number of analyzed observations are given for each star. 

Periods longer than 10 days were found only for 11 stars, i.e., about 30%, of all 
detected periodic variables. We believe that such a low proportion of slowly rotating 
stars is a result of observational selection: the expected amplitudes of periodic varia- 
tions of high activity stars are substantially higher, hence easier to detect than in case 
of inactive stars. Several more M type stars with short rotation periods of the order of 
one or a few days are known in the literature (Messina et al. 2003) but we found only 
two stars with long rotation periods: GJ 411 with /Vot = 48 d (Noyes et al. 1984), and 
GJ 699 (Barnard star) with Prot » 130 d (Benedict et al. 1998). We include both stars 
into the further analysis. 

Notes on individual stars. 

GJ 84 is a visual binary with the secondary component lying 0''44 from the primary 
(Golimowski et al. 2004). The luminosity difference is 4.6 mag in the Fl lOW filter at 
the NICMOS camera, so the contribution of the secondary is negligible, both in the 
V-band and (most likely) in X-rays. Aliases 0.9749 d and 1.020 d are rather excluded 
due to the low X-ray activity of the star. 

GJ 103 (CC Eri) is a spectroscopic binary with components of spectral types K7 + 
M3 and the combined type MO (NEXXUS). The photometric period is the same as the 
orbital period (Evans 1959). 

GJ 1054A is a component of double line binary with an unknown period (Gizis, 
Reid and Hawley 2002) which, however, according to the authors, should be short. 
The most significant period in the ASAS photometry is given in Table 1 (0.513 d, 
©AoV = 15.4), its aUas 18.86 d is also significant, (0aoV = 13.3) but we consider it as 
less likely due to very high activity of the star. 

GJ 182 is a single flare star with a rotation period known from literature of 4.565 d 
(Byrne et al. 1984). We have not found any significant variability at this period. Max- 
imum value of ©AoV = 26.4 is at a period of 1.2921 d for the 106 data points between 
hjd 2496.9 and 31 10.5, however we consider the period of 4.410 d with ©aoV = 19.9, 
obtained for the whole data set, as the most probable. 

GJ 2036A was observed together with GJ 2036B by ASAS. Component B is at an- 
gular distance 8"and is fainter by about 1 mag in the V-band. The total V magnitude of 
the binary as measured by ASAS is 10.75 mag. We attribute the photometric variabiUty 
to the brighter component. 

GJ 3331 A is the brightest component of the triple system. The pair BC is 8''2 
from the primary and the differences in magnitudes between components A and BC 
are AV = 0.722, AR = 0.574, A/ = 0.220 (Jao et al. 2003). We assume that photometric 
variability is related to GJ 3331 A. Its parallax measurement has a large error (NEXXUS 
database) so we do not include this star in Table 3. 

GJ 205 - aUases at 0.971 d and 1.028 d are rather excluded due to low activity of 
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9.00 
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Fig. 1. Photometric measurements of periodic stars vs. hjd = HJD -2 450 000 d (left panels). Observations 
showing periodic variations are plotted as filled symbols. Right panels show filled symbols folded with the 
Usted period. The brightness at right panels is corrected for seasonal variations of luminosity. The number 
of filled symbols and the total number of analyzed observations are also given. 
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the star. 

GJ 375 is a double line spectroscopic binary with a period of 1.875 d and nearly 
identical components (Montes et al. 2006). Photometric period is identical with the 
spectroscopic one. Mass, given in Table 1, refers to one component. 

TWA 5 - appears to be a multiple system. This is a young system belonging to 
the TW Hya association. The main component consists of a spectroscopic binary and 
the third component is separated by 0"054 (Brandeker et al. 2003). All three stars are 
of spectral type Ml. 5. In addition, there is also a brown dwarf at angular distance 
of 2". Based on XMM observations Argiroffi et al (2005) give log(Z,x) = 29.8 and 
log [Lx /L\io\ ) = — 3 . 1 . for the triple M dwarf system. The photometric period found in 
the ASAS data may probably be related to the orbital period of a spectroscopic binary. 
Another possible period is an alias 3.512 d. The mass given in Table 1 is based on the 
absolute magnitude of the total system and, therefore, very likely overestimated. 

GJ 494 (DT Vir) has a low luminosity companion at angular separation of 0"48 
(February 2000, Beuzit et al. 2004), fainter in the K band by 4.4 mag, which makes its 
contribution completely negligible, both in V and in X-rays. 

GJ 55 1 = Proxima has the lowest luminosity in our sample of stars with detected 
periods. ASAS photometric data indicate a clear periodicity of 82.5 days, very close to 
the rotation period found by Benedict et al. (1998). We have not found any significant 
signal at a period of 42 days, visible in part of the data analyzed by Benedict et al. 
(1998), or around 30 days, reported by Guinan and Morgan (1996). 

GJ 618A has a nearby companion (5 "4) of spectral type M5, which is fainter in 
the y-band by 3.5 mag. This makes its photometric contribution completely negligi- 
ble and we can safely assume that the photometric variability is related to the brighter 
component. Nevertheless, the X-ray luminosity of M5 type star may contribute signif- 
icantly to the total X-ray emission of the binary if the star is very active, so the value 
of Rx = 4.7, listed in Table 1, should be treated as an upper limit for the X-ray flux of 
GJ618A. 

GJ 2123A is a very active M3 star with an M4 companion about 1 .8 mag fainter in 
the V-band, lying at angular distance of 2"8. We attribute the photometric variability 
to the brighter component. 

GJ 669 A has a companion at angular distance of 16", fainter by about 1.6 mag in 
the V-band. Both stars are unresolved in the ROSAT and ASAS data. 

GJ 674, there are also significant peaks of AoV statistics at aliased periods of 
0.9709 d and 1.028 d with ©aov « 10. We prefer 33.29 d due to low activity of the 
star. 

GJ 799 (AT Mic) is a visual binary with angular separation between components 
of 3"3, too close to be resolved either by ROSAT or by ASAS. Both components have 
almost the same luminosity in V-band and we assume that their X-ray luminosity is 
also the same. Values listed in Table 1 are related to a single component. 

GJ 803 (AU Mic) is a young, single star with rotation period of 4.865 d (Torres et 
al. 1972). We see the strongest signal at a period of 4.8482 d in the ASAS data (62 data 
points between hjd 3450.9 and 3718.5), but there is also a significant power at 4.8604 
and 4.8617 d for the whole data set. This may be related to the differential rotation of 
the star. 

GJ 1264A has a companion at a distance of 1''3, fainter by 1.2 mag in the V-band. 
The stars are unresolved in ROSAT and ASAS data. We assume that all activity and 
variability is related to the brighter component. 

GJ 841 A is a spectroscopic binary with a period of 1.1248 d (Jeffries and Bro- 
mage 1993). A similar period of 1.1237 d is visible in the whole ASAS data set (with 
the AoV statistics equal to 12.4). More significant (aliased) period of 0.5283 d with 
©AoV = 12.5 does not fit to the Jeffries and Bromage data. 
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GJ 867A (FK Aqr) has a companion separated by 23'! 5 and fainter by 1.9 mag 
in the V-band. The stars are unresolved in ROSAT and ASAS data. GJ 867A is a 
spectroscopic binary with an orbital period of 4.083 d (Herbig and Moorhead 1965). 
Bopp and Espenak (1977) also found a photometric period of 4.08 d, but Byrne et 
al. (1987) claim a photometric period of 4.39 d. Cited above photometric periods are 
close to, but not identical with our value 4.232 d. The mass given in Table 1 is likely 
overestimated. 

GJ 890 (HK Aqr). Barnes and Collier Cameron (2001) give a rotation period of 
0.4307 d for this star. We obtained a similar period of 0.4311 d. 

GJ 897A has a companion separated by 0"5 and fainter by 0.5 mag in the V-band. 
The stars are not resolved by ROSAT and ASAS. We attribute the photometric variabil- 
ity to the brighter component. 

3.2 Empirical 'Rirnover Times 

To derive the empirical turnover times for low mass stars the procedure applied by 
Stgpiefi (1994, 2003) was used, except that the analyzed stars were grouped in the 
present paper according to mass rather than B — V. Natural breaks in the mass distribu- 
tion were used for this grouping. Two stars: GJ 182 and TWA 5 A stand out as the most 
massive stars. Because no slowly rotating stars with similar masses occur in Table 1, 
the stars were excluded from the further analysis. 10 "high mass" stars have masses be- 
tween 0.56 M© and 0.65 M©, 10 "medium mass" stars have masses between 0.43 M© 
and 0.52 M© and 8 "low mass" stars have masses between 0.33 M© and 0.39 Mq. 
Two stars with masses 0.2 M© and one with 0.11 M© are termed "very low mass" stars 
(Table 2). 



Table2 

Values of the coefficients a and b from Eq. (4) for four mass intervals and the resulting turnover times 



stars 


av. mass 


No. of stars 


a 


b 




high mass 


0.61 


10 


-2.94 ±0.08 


-0.050 ±0.006 


30±4 


medium mass 


0.47 


10 


-3.02±0.18 


-0.040 ±0.007 


38±6 


low mass 


0.37 


8 


-3.14±0.08 


-0.029 ±0.003 


53±6 


very low mass 


0.17 


3 


-3.00 


-0.016 ±0.005 


95 ±30 



Fig. 5 presents values of \ogRx vs. /rot for all four groups of stars. Different groups 
are marked with different symbols. Linear relations of the form 

\ogRx = a + bP,oi (4) 

are fitted to the data within each group. For the first three groups both coefficients were 
determined with the least squares method. Because only one short period star with very 
low mass is plotted in Figs. 1^ the value of a in this case was assumed to be equal 
to —3. Table 2 gives values of the coefficients from Eq. (1) for each group, together 
with their errors. The value of the coefficient b for very low mass stars is extremely 
uncertain due to a low number of data. 

Values of empirical turnover time can be found from relation: =—cb^^, where the 
normalization factor c remains undetermined unless one wants to compare the obtained 
values of lc to values obtained with another method (for details, see Stgpien 1994). A 
value of the normalization coefficient c = 1.515 was found by fitting the presently 
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150 



Fig. 5. X-ray to bolometric flux ratio vs. rotation period for stars from Table 1. Straiglit lines are the least 
square fits to observational data grouped according to mass. Asterisks correspond to high mass stars, pluses 
to medium mass, diamonds to low mass and triangles to very low mass stars. Average mass for each group 
is indicated at the fitted line. 



calculated turnover times, supplemented with the values obtained by Stgpien (1994) 
from the X-ray data of G-K type stars (B — V between 0.55 mag and 1 .25 mag), to 
the relation determined by Stgpien (2003) from Call emission. The latter relation was 
found by using correctly reduced values of the net calcium emission flux, contrary 
to, so called, excess calcium emission flux used earlier by St^pien (1994). The last 
colunm of Table 2 gives normalized values of Xc and Fig. 6 shows Xc v*. mass. The 
relation from Stgpieri (2003) is plotted as a solid line and the values of turnover times 
found from X-ray data are given with their errors. The older data are plotted with 
B — V values transformed to stellar mass. It is seen from Fig. 6 that, after passing a 
plateau for masses in the range 0.8 M0-O.6 M©, turnover time increases again for still 
lower masses. The increase is not as steep as resulting from the extrapolation of the 
theoretical values (Kim and Demarque 1996). Overplotted are values of turnover times 
found by Pizzolato et al. (2003). They used a different sample of stars and a somewhat 
different method but the agreement between both sets of data is satisfactory. They 
were not able to determine turnover times for stars less massive than 0.6 M© due to a 
shortage of data for slowly rotating stars. 

Fig. 7 shows the X-ray data of the stars from Fig. 5 plotted vs. Rossby number 
which was calculated for each star using turnover time from Table 2. As expected, 
the scatter is considerably lower than in Fig. 5 and the exponential relation between 
Rx and Ro seems to describe satisfactorily the fit \og{Rx) = —(3.05 ±0.07) — (1.45 ± 
0.12)Ro. Similar relation has been shown to hold for G and K type stars (Stgpien 1994: 
log{Rx) = -(3.71 ±0.27) - (1.46±0.13)/?o). The free term differs somewhat between 
both relations but the slope agrees very weU. 

3.3 Age-Period and Period-Activity Relations 

Active stars from the lower MS are expected to lose angular momentum via a mag- 
netized wind. Observations of stellar clusters show indeed that single stars spin down 
with age. Unfortunately, accurate measurements of rotation periods are available for 
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Fig. 6. Empirical turnover times with error bars, determined from X-ray fluxes, vs. stellar mass. Data from 
Table 2 are supplemented with older determinations for more massive stars (St?pien 1994). Values of the 
turnover times are scaled to data found from calcium emission fluxes (soUd Une, Stgpieii 2003). Diamond 
give values of empirical turnover times from Pizzolato et al. (2003). 



members of young clusters only. In case of M 67, the oldest cluster (4 Gyr) for which 
useful measurements exist, no more than just vsini values have been measured for 
some stars. The available data may be sufficient to calibrate the age-rotation relation 
for solar type stars but M type stars evolve substantially slower and we need accurate 
data for stars much older than members of the investigated clusters. Such stars are ex- 
pected to exist among field stars in the solar vicinity but accurate determination of field 
star age is difficult. We are forced to use proxies only statistically related to age. A 
commonly used proxy of this kind is kinematic class (Eggen 1969) based on the (U, V) 
plane, where U, V and W are the components of space motion of a star. Table 3 gives 
star names, components of proper motions in right ascension and declination, radial 
velocity, the calculated values of the U, V and W components for stars from Table 1 for 
which complete data, including radial velocity, are available. Proper motions and radial 
velocities are from the NEXXUS database, U-axis is directed opposite to the Galactic 
center, V-axis - in direction of Galactic rotation, W-axis - in direction of the North 
Galactic Pole. YD means Young Disk with kinematical parameters —20 <U < 50, 
30 < y < 10, -25<W< 10, OD - Old Disc, YD/OD Young or Old Disc, OD/H - 
Old Disk or Halo. 

Known spectroscopic binaries with short periods were excluded because their ac- 
tivity levels and variability periods (assuming synchronization of orbital and rotational 
period) are not related to ages. The last column gives kinematic class. We identify 20 
stars belonging to young disk (YD), 1 star falling at the border between young and old 
disk (YD/OD) and 4 stars belonging to old disk (OD). Both stars added to our sample 
from Uterature belong to old disk (Fleming et al. 1995), so we have all in all 6 OD 
stars, although GJ 699 may even belong to halo (Delfosse et al. 1998). 

As it is seen from Table 3, all short period stars belong to YD or YD/OD. However, 
there are also 4 YD stars in our sample with periods longer than 10 days, including 
Proxima with a period of 82.5 d - the longest one found in our search. Such a mixture 
of rapidly and slowly rotating stars is not surprising. YD stars are on average about 
3-4 Gyr old (Meusinger, Stecklum and Reimann 1991) but individual stars may have 
ages up to several Gyr. If Proxima has a connmon origin with a Cen (Wertheimer and 
Laughlin 2006), its age is about 6 Gyr, yet its kinematical characteristics are typical of 
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young disk. 

The average period of YD stars is 9.4 d and the median value is about 2 d. Contrary 
to YD stars, we do not have any rapidly rotating stars among OD objects. All six stars 
have periods longer than 30 d with a mean value of 58.8 d and a median of 47 d. We 
assume that OD stars have age of the order of 10 Gyr (Meusinger et al. 1991). 

Comparison of both median values indicates that angular momentum of M type 
stars decreases, on average, by a factor of 25 in 10 Gyr but this factor may vary between 
5 and 100 if the initial rotation periods of individual stars are taken from observations 
of young clusters (Barnes 2003). 

Based on observations of csin/ of solar type stars Skumanich (1972) found an em- 
pirical relation describing the dependence of the rotation period on time: P^i °= t^^^- 
More numerous observations of stars in several young clusters showed that the Sku- 
manich law breaks down for rapidly rotating stars which show the effect of saturation, 
visible as a flattening of Pmt{t) when t ^Q. Considering slowly rotating stars we still 
may ask, how universal the Skumanich law is? Is it mass independent? A number of 
formulas describing spin down rate of single stars have been suggested in the litera- 
ture. One of them, derived by Kawaler (1988), has a form: dco/df = Aco^, where the 
coefficient A depends on stellar mass and radius. Neglecting time variability of stellar 
parameters we recover the Skumanich law: coocAf Here, the dependence on stellar 
parameters appears only through the coefficient of proportionality A whereas the func- 
tional dependence on time {t~^l^) is mass independent. Sills et al. (2000) assumed that 
the Kawaler formula can be applied even to VLMS. Another spin down formula, in 
which mass dependence occurs only via an exponential term containing turnover time, 
was derived by St§pieii (1988, 2006) 

-^=7xlO-W^«/«-335 (5) 
at 

where CO is expressed in 1/d, and t in years. After integration the above formula predicts 
different functional dependence of the rotation period on time for different turnover 
times, hence masses. In particular, stars with turnover times substantially longer than 
the Sun spin down faster and rotate slower at the same age. Fig. 8 shows the time vari- 
abiUty of a rotation period for stars with turnover times from Table 2 and (assumed) 
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Tables 
Kinematic properties of stars 











u 


V 


W 


remarks 




[nifls] 


[nifls] 




[km/s] 


[km/s] 


[km/s] 




GJ84 


1317.5 


-173.9 


23.3 


45.2 


-44.5 


-6.8 


OD 


GJ 176 


659.8 


-1114.7 


26.0 


22.6 


-57.4 


-14.7 


OD 


GJ 182 


37.2 


-93.9 


18.2 


11.3 


-16.8 


-9.2 


YD 


GJ 2036A 


132.9 


73.9 


40.6 


8.4 


-35.1 


-20.3 


YD/OD 


GJ205 


763.0 


-2092.8 


8.5 


-22.2 


-55.6 


-10.2 


OD 


GJ 3331A 


49.1 


-32.7 


20.0 


11.5 


-14.1 


-9.0 


YD 


GJ 3367 


-85.2 


-77.3 


-23.0 


-24.4 


8.0 


-6.8 


YD 


GJ358 


-526.1 


356.9 


7.0 


28.5 


-6.7 


-2.8 


YD 


GJ382 


-152.9 


-242.9 


7.6 


2.4 


-12.5 


-2.8 


YD 


GJ431 


-715.8 


170.8 


6.0 


32.8 


-17.2 


-1.2 


YD 


GJ494 


-618.8 


-16.6 


-11.7 


29.3 


-17.2 


-10.3 


YD 


GJ551 


-3775.4 


769.3 


-21.6 


28.9 


1.5 


13.7 


YD 


GJ 569A 


276.1 


-122.4 


-7.2 


-7.8 


3.2 


-13.3 


YD 


GJ 9520 


78.4 


129.4 


6.5 


-0.4 


9.5 


4.3 


YD 


GJ618A 


-729.3 


991.2 


43.0 


-35.6 


-3.2 


54.9 


OD 


GJ 2123A 


65.8 


-109.4 


-18.0 


17.8 


1.3 


-9.0 


YD 


GJ 669 A 


-221.2 


349.0 


-33.6 


34.5 


-19.1 


-3.6 


YD 


GJ674 


573.4 


-879.6 


-10.0 


14.5 


-5.1 


-19.3 


YD 


GJ729 


637.6 


-192.5 


-11.5 


13.0 


-1.2 


-7.1 


YD 


GJ 799 AB 


269.3 


-365.7 


5.0 


2.7 


-15.5 


-16.2 


YD 


GJ 803 


280.4 


-360.1 


-4.6 


10.2 


-16.4 


-10.4 


YD 


GJ 1264A 


300.8 


-291.1 


-4.5 


26.0 


-19.2 


-1.7 


YD 


GJ 867A 


450.6 


-79.9 


-8.7 


17.8 


-10.3 


-2.1 


YD 


GJ890 


101.6 


-24.1 


7.0 


6.2 


-3.1 


-10.9 


YD 


GJ 897A 


317.0 


-234.0 


-2.6 


13.1 


-24.6 


-7.8 


YD 


GJ411 


-580.5 


-4770.0 


-84.7 


-46.2 


-53.8 


-74.3 


OD 


GJ699 


-798.7 


10337.8 


-110.6 


141.1 


4.5 


18.2 


OD/H 



The columns give name of a star, proper motion in right ascension and declination, radial 
velocity, three galactic components of velocity, kinematical class (see text) 



initial rotation period of 1 d. For comparison, Pmi{t) of a star with turnover time cor- 
responding to the solar value (15.5 d) and to the value from plateau (25.5 d) are also 
given. Overplotted are the observed rotation periods of the Sun, Proxima and all six 
OD stars (plotted at an age of 9.9 Gyr for better visibility). The observational point at 
an age of 4 Gyr corresponds to a period of 17 d, calculated from the average value of 
V sin i of four solar type starswith0.60mag<fi — y<0.69mag measured by Pace and 
Pasquini (2004). Two highest lying curves correspond to Xc = 95 d and to this turnover 
time minus its estimated error, /.e., 1^ = 65 d. It seems that the lower value better agrees 
with observations of the least massive stars. The overall agreement with observational 
data supports the slow down model described by Eq. (5). Note that the rotation period 
of a star with the initial period of a fraction of a day would vary along the curve lying 
correspondingly lower in Fig. 8 than the curve starting at Pim = 1 d. 

Extensive observations of chromospheric emission Unes and X-ray fluxes of so- 
lar type stars in many open clusters resulted in well determined age-activity relations 
(e.g., Barry 1988, Soderblom, Duncan and Johnson 1991, Haisch and Schmitt 1996, 
Giidel, Guinan and Skinner 1997). However, a direct extrapolation of a relation cal- 
ibrated on stars with masses 0.7-1.1 M,:, to M stars with masses several times lower 
is not obvious a priori and needs to be examined. Observations of M stars of a given 
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Fig. 8. Time evolution of stellar rotation period, resulting from integration of Eq. (5) with the initial value 
of 1 d. Different curves correspond to different turnover times, as indicated. Overclouded are observational 
data with the following coding: diamond - the average rotation period of four solar type members of M 67, 
calculated from their v sin values, open circle - the Sun, cross - Proxima and asterisks - six old disk stars 
from Table 3. 



age show that rotation rate and activity level increase towards later spectral subtypes 
(Delfosse et al. 1998, Scholz and Eisloffel 2005) but a mean level of activity decreases 
with age for a given spectral type (Hawley, Reid and Tourtellot 2000). Similarly as in 
case of rotation data, the activity measurements of stars with well determined ages are 
restricted to nearby open clusters and the observations of activity level of still older M 
type stars with ages of the order of 10 Gyr refer to field stars. Two such searches have 
recently been performed: one in Ha (Silvestri, Hawley and Oswalt 2005) and the other 
in X-rays (Feigelson et al. 2004). Silvestri et al. (2005) observed 139 M stars in wide 
binaries with white dwarf companions. Ages of the binaries were determined from the 
cooling age of the white dwarf components. Most of the investigated M dwarfs have 
ages below 4 Gyr but there are also many stars with ages up to 10 Gyr (regretfully the 
authors do not list ages of individual stars). The results confirm earlier findings that the 
mean activity level increases with decreasing mass. The analysis of kinematics of the 
investigated stars showed that active stars (of dMe type) have lower space velocities and 
are confined to thin disk. Inactive stars (of dM type) have systematically higher space 
velocities, characteristic of thick disk. Feigelson et al. (2004) analyzed the results of 
a deep, pencil-beam X-ray observations in a Galactic high-latitude field, obtained with 
Chandra. The observations of late-type MS stars are compared with predictions based 
on convolution of X-ray luminosity function with the known spatial distribution of old 
disk stars. The authors came to the conclusion that stellar X-ray luminosities decrease 
with age like Lx °^ t~^ over 1 < / < 11 Gyr, rather than Lx°^t^^ resulting from simple 
power laws: Prot °^ t~^l'^ (Skumanich 1972) and o= Ujot (Pallavicini et al. 1981, Piz- 
zolato et al. 2003). The result shows that at least one the above power laws does not 
apply to M stars which dominate stellar source counts in the search by Feigelson et al. 
(2004). 

Stgpieri (1988, 1994) suggested an alternative period-activity relation, in which the 
activity level depends exponentially on the Rossby number. 

We have shown in the previous section that the X-ray activity of M stars decreases 
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with the Rossby number similarly to that of G-K type stars i.e., 



logRxO^ 



-(1.45±0.12)/?o. 



(6) 



Quality of fit of both, power and exponential relations to observational data is sim- 
ilar (Hempelmann et al. 1995) so it is not possible to discriminate between them on the 
basis of the rotation-activity diagram alone. However, the age-activity relation result- 
ing from Eq. (6), together with Eq. (4), is in a much better agreement with the results 
obtained by Feigelson et al. (2004) than the power relation Lx°^t^K 

Using turnover times given in Table 2, Rossby numbers of six old disk stars were 
calculated and the average value obtained: /?o=1.2±0.4. The uncertainty of this value 
comes from errors of turnover times and the scatter of individual Rossby numbers. 
Assuming that the obtained value of Ro is typical for old M type stars and that young 
M type stars have /?o -C 1 , we obtain: 



X-ray observations of Hyades (Reid et al. 1995) and Praesepe (Franciosini et al. 
2003), open clusters with ages about 600 Myr and 800 Myr respectively, show their M 
stars still in saturation regime. Assuming the age of M stars at the end of saturation 
regime « 1 Gyr and of old stars « 10 Gyr, the 10-fold increase in age should result, 
according to Eq. (7), in a decrease of X-ray flux by a factor 10~' Z This prediction 
can directly be compared with observations. All six OD stars were observed in X-rays 
(see Table 1) and their average value of logRx is equal to —4.84. The average value 
of logRx of young, rapidly rotating stars is about —(3.0—3.2) (Pizzolato et al. 2003, 
see also Fig. 5), so the observed difference (log^x)young — (log^x)oid ~ 1.6—1.8 is in 
excellent agreement with the predictions based on Eq. (7) and in a fair agreement with 
the results by Feigelson et al. (2004). 

Based on the results of the present paper we can explain why M type stars show 
higher activity level even if they are older and rotate much slower than more massive 
stars. As it was shown, the main parameter controlling activity is the Rossby number, 
not rotation period. The results of the present paper show that the empirical turnover 
time increases with decreasing mass, at least for early M types. Rotation periods of 
the old disk M type stars are not much longer than their turnover times so the resulting 
Rossby numbers are of the order of unity. Such stars should be considerably more 
active than massive stars with the same rotation periods but shorter turnover times. 
Our analysis applies to early M type stars up to about M4, with masses above the 
full convection Umit. There is indication that empirical turnover time increases further 
for stars with masses as low as 0.17 M© but a significantly more numerous sample is 
needed to draw meaningful conclusions. 

4 Conclusions 

We checked 180 nearby X-ray active M type stars for periodic light variations. Anal- 
ysis of 6 seasons of ASAS photometric observations led to the detection of periodic 
variations in 31 objects with typical half-amplitude of 10 mmag and periods between 
a fraction of a day and 83 d. We interpret these variations as resulting from rotational 
modulation of non-uniformly spotted stellar surface. Prior to our search only three M 
type stars with long (at least a few tens of days) rotation periods were known. We 
detected 9 more stars with rotation periods longer than 10 d which is the approximate 
limiting period of saturation for early M dwarf stars. Our data show that also M type 
stars obey period-activity relation such that activity decreases with increasing period. 



(l0g/?x) 



young 



(log/?;,)oid = 1.7 ±0.50. 



(7) 
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The relation is mass dependent, similarly as in G and K type stars. To obtain a single 
period-activity relation rotation for all stars, their periods should be scaled down by 
mass dependent factors (turnover times). The analyzed stars were divided into four 
mass intervals and turnover times were determined for three of them. The fourth bin, 
containing the least massive stars has too few data points to obtain its reliable value 
but a continuous increase of the turnover time with decreasing mass, down to masses 
of 0.1-0.2 Mg,, is suggested. We demonstrated that the formula for spin down rate 
of single stars, derived by Stgpiefi (1988) also correctly predicts the evolution of ro- 
tation period of M type stars when newly determined turnover times are used. The 
analysis of space motions showed that all OD stars have long rotation periods with a 
median value of 47 days. This clearly demonstrates that M type stars spin down as 
they age. Assuming 10 Gyr for the age of OD stars we showed that the exponential 
period-activity relation, together with the spin down rate formula, correctly predicts 
the decrease of X-ray activity with age, as observed by Feigelson et al. (2004). The 
observations of Feigelson et al. are at odds with an age-activity power relation result- 
ing from the Skumanich law and period-activity relation suggested by Pallavicini et al. 
(1981) and Fizzolato et al. (2003). 

It would be desirable to extend the period search to late M stars but to this purpose 
photometric observations reaching fainter stars than present day ASAS are needed. 
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